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In the present report, let us adress the issues related to simulation of decays for particle embodied in full 
production and decay chains of Monte Carlo programs set-up for experiments such as at LHC or LC. Both 
technical issues related to the way how the events may be stored in event records and issues related to physics 
(in particular non-factorizable correlations of the Einstein-Rosen-Podolsky type) will be reviewed on the basis of 
practical examples. We will limit our discussion to the case of r lepton and W boson decays, but similar problems 
(and solutions) may arise also in case of simulation for other intermediate states or particles. Examples related 
to construction of physics observables will be also given. In particular the method of measuring the CP parity 
properties of the h — tt coupling at LC will be explained. 
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1. Introduction 

Since many years, intensive studies are being 
performed to design future software architectures 
for experiments on proton proton colliders, such 
as the Tevatron or the LHC ^ and high energy 
e+e" linear colliders such as JLC, NLC T or 
TESLA [g. 

One of the important ingredients in such de- 
signs is the data structure for storing the Monte 
Carlo events. It is generally accepted that the 
data structures based on objects such as parti- 
cles, clusters, strings, etc. with properties such 
as tracks, momenta, colour, spin, mass, etc. and 
on the relations explaining the origins and descen- 
dants of the objects is the most convenient one. 
This is the case at present j4j , and it is also envis- 
aged for the future, see [S]. At the same time such 
a picture is in conflict with the basic principles of 
quantum mechanics. Einstein-Rosen-Podolsky 
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paradox is an example of such phenomena. A 
general problem is that the quantum state of a 
multiparticle system cannot (at least in principle) 
be represented as a statistical combination of the 
states defined by the products of the pure quan- 
tum states of the individual particles. It is thus 
of the utmost importance to examine whether the 
approximation enforced by the data structure is 
purely academic, or if it rather represents a real 
difficulty, which may affect the interpretation of 
the future data. In fact in some cases alternative 
methods can be designed and are in fact used as 
well. 

It would not be a serious problem if the pre- 
dictions of the Standard Model used in the in- 
terpretation of the future data could be provided 
by a single program, black box, without any need 
of analysing its parts. Then anything that would 
be measured beyond the prediction of such a hy- 
pothetical Monte Carlo program would be inter- 
preted as "new physics" . Agreement, on the other 
hand, would constitute confirmation of the Stan- 
dard Model, as it is understood at present (and 
proper functioning of the detector as well). How- 
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ever, even in such an extreme case it is very useful, 
for the purpose of experimental studies, to manip- 
ulate with the terms responsable for the signature 
of the 'new physics'. In this way experimental 
strategies can be refined, if the new effects can 
be placed in well defined and phenomenologically 
simple modules. 

Because of the complexity of the problem, 
Monte Carlo predictions need to be dealt with 
by programs describing: the action of the detec- 
tor and of the analysis, on the experimental side, 
and various effects, such as those from hard pro- 
cesses, hadronization, decay of resonances, etc., 
on the theoretical side. Every part is inevitably 
calculated with some approximation, which need 
to be controlled. 

2. Event record and decay interface 

In the first part I will discuss solution we 
used in KORALZ |2| - the program widely used 
at LEP for the simulation of r-lepton pair pro- 
duction and decay, including spin and QED 
bremssstrahlung effects. Even though spin effects 
are non-treatable in the scheme where properties 
are attributed to individual particles only, it is 
the very method used there. As described in ref. 
[S] the algorithm of spin generation for any indi- 
vidual event was consisting of the following steps: 

1. An event consisting of a pair of r leptons, 
bremsstrahlung photons, etc., was gener- 
ated. 

2. Helicity states for both t+ and t~ were 
generated. At this point, an approximation 
with respect to quantum mechanise was in- 
troduced. 

3. Information on these helicty states, includ- 
ing the definition of quantization frames, 
i.e. the relation between t's rest frame and 
laboratory frame, was then transmitted to 
TAUOLA |7I8I9| . the package for the genera- 
tion of r-lepton decays. 

4. Finally TAUOLA performed decays of 100% 
polarized r's, and the event in the HEPEVT 
common block was completed. 



The solution for the spin treatment of r leptons 
at LEP was optimal. On one side, a convenient 
picture of particles with properties, origins and 
descendants could be used and, on the other, a 
complete full spin solution |10lll| was available, 
if necessary. 

Let us now turn to another example of the 
spin implementation algorithm. It is taken 
from ref.^21- The algorithm, essentially that of 
KORALZ, was adopted to work with any Monte 
Carlo program providing the production of r- 
Icptons. If the generated events are stored in the 
format of a HEPEVT common block, then the algo- 
rithm consisting of the following basic steps can 
be used: 

1. Search for r-leptons in a HEPEVT common 
block (filled by any MC program). 

2. Check what the origin of r-lepton is: 
Z, 7, W, h, or eventually, 2 — > 2-body 
process such as: e^e~, {uu), {dd)^T^T^. 

3. For the 2 2-body process of r-pair pro- 
duction, it is sometimes possible to calcu- 
late the r polarization as a function of the 
invariant mass of the r-lepton pair and an- 
gle between the directions of r-leptons and 
incoming effective beams (in the rest frame 
of r-pair). 

4. If in addition to the r-leptons, photons or 
partons (gluons, quarks, etc.) are stored in 
HEPEVT common block, one needs to define 
the "effective incoming beams" . 

5. From such an information one can generate 
r helicity states and define the relation be- 
tween the r rest frame and the laboratory 
frame. Optionally complete spin effects can 
be implemented as well, see |18| . 

6. The r decay is generated with the help of 
TAUOLA and HEPEVT common block is ap- 
pended with the r's decay products. 

Leading spin effects are nicely reproduced by 
the above set of programs. A more complete dis- 
cussion can be found in ref. |12j . 

Let us stress, that the presented above solu- 
tion, require certain minimal discipline in a way 
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Figure 1. Typical relations in present day event 
record. Daughter index JDAHEP points upward in 
a tree, however mother index JMDHEP may point 
to the previous copy of the event as well. The sta- 
tus code ISTHEP is ocasionally used like a pointer 
toward higher level of the event record listing. As 
a consequence in the event record there may ap- 
pear ambiguities and/or loops for searching algo- 
rithms. 



how the event records are filled in. At present, 
there is a strong tendency to store in event record, 
not only information on the 'real particles', but 
also on the results of the simulation at the par- 
ton shower level as well as of the hard process 
alone. As a consequence not only the same en- 
tries for the otherwise well defined partilces like 
r-leptons are duplicated or even triplicated, but 
also the relation between different entries is not 
anymore of the tree-type and links between par- 
ticles are not reversible. The link upwards from 
A to B docs not mean that there is downward 



link from B to A. This creates multitude of trou- 
bles for the algorithms analyzing events, see fig.ffl 
However at present our solution for interfacing 
decay packages with the host programs based on 
HEPEVT event record as filled by PYTHIA and 
HERWIG seem to work in all cases studied by us 
|TT^) . For the case of complete spin correlations, 
we found it more convenient to abandon direct 
use of spin information provided by the host pro- 
grams. Instead we choose to calculate complete 
density matrix anew, from the kinematical con- 
figuration provided by the host program. 

3. Higgs boson parity measurement 

Let us sketch the basic principles behind the 
proposed measurement, in the case, when simul- 
taneosly scalar and pseudoscalar couplings are al- 
lowed in Htt vertex 

/i fiV(cos0 -I- «sin(/)75)T. (1) 

If non-zero CP-odd admixture to the Higgs is 
present, the distribution of the Higgs production 
angle is modified I14I15I16] . We have simulated 
production angular distributions as in the SM, 
but this assumption has no influence on the va- 
lidity of the analysis. In order to study the sen- 
sitivity of /i — > T+T^ observables, we assume a 
SM production rate inependent of the size of the 
CP-odd admixture. 

The production process e^e~ Zh —>■ 

fj,^ {qq)T^ T~ has been chosen, as an represen- 
tative example, and simulated with the Monte 
Carlo program PYTHIA 6.1 f7\ The Higgs bo- 
son mass of 120 GeV and a centre-of-mass energy 
of 350 GeV was chosen. The effects of initial state 
bremsstrahlung were included. For the sake of 
our discussion and in all of our samples the r de- 
cays have been generated with the TAUOLA Monte 
Carlo library 9 8 7,. As usual, to facilitate the 
interpretation of the results, bremsstrahlung ef- 
fects in decays were not taken into account. Any- 
way, with the help of additional simulation, we 
have found this effect to be rather small. To 
include the full spin effects in the h — » r+r^, 
P^Dt{i't)) TT^TT^ decay chain, the 

interface explained in Ref. ^Sl was used. 
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The Higgs boson parity information must be 
extracted from the correlations between r"*" and 
T" spin components, which are further reflected 
in correlations between the r decay products in 
the plane transverse to the r+r" axes (19ii20) . To 
better visualize the effect, let us write the decay 
probability, using the conventions of Ref. : 

r(;i™x ^ T+T-) ^ +si^i?(20) si, [2) 

where R(2(j)) can be understood as an operator 
for the rotation by an angle 2(j) around the || di- 
rection. The s'^ and s'^^ are the polariza- 
tion vectors, which are defined in their respective 
rest frames. The symbols ||/_L denote components 
parallel/transverse to the Higgs boson momen- 
tum as seen from the respective rest frames. 

The method relies on measuring the acopla- 
narity angle of the two planes, spanned on 
decay products and defined in the p'^ p~ pair rest 
frame. The acoplanarity angle (/?*, between the 
planes of the and p~ decay products is defined. 
The angle is defined first, with the help of its co- 
sine and two vectors n-t normal to the planes, 
namely n± = p^± x p^o, costp* = |"+||"^| . 

To distinguish between the two cases: ip* and 
27r — ip* it is sufficient, for example, to find the 
sign of p^- -11+. When it is negative, the angle ip* 
as defined above (and in the range < (p* < tt) 
is used. Otherwise it is replaced by 2tt — tp* . 

Additional selection cuts need to be applied. 
The events need to be divided into two classes, 
depending on the sign of yij/2, where 

E^+ + £"^0 -Ett- + E^^o 

The energies of tt^ , tt" are to be taken in the re- 
spective rest frames. In Refs. |19I20) the meth- 
ods of reconstruction of the replacement t"^ rest 
frames were proposed with and without the help 
of the r impact parameter. We will use these 
methods here as well, without any modification. 

To test the feasibility of the measurement, some 
assumptions about the detector effects had to be 
made, see refs. |19I2(J| for more details. 

3.1. Numerical results 

We have used the scalar-pseudoscalar mixing 
angle (j) = ^. In Fig. |21 the acoplanarity distribu- 



tion angle (p* of the p^ p^ decay products which 
was defined in the rest frame of the reconstructed 
p^ p^ pair, is shown. The two plots represent 
events selected by the differences of tt^tt" ener- 
gies, defined in their respective rest frames. In 
the left plot, it is required that yiy2 > 0, whereas 
in the right one, events with yiy2 < are taken. 
This figure quantifies the size of the parity effect. 
The size of the effect is substantially diminished 
when a detector-like set-up was included for 
rest frames reconstruction, in exactly the same 
proportion as in Ref. |12j, nonetheless parity ef- 
fect remain visible. 

The fitting procedure was repeated 400 times 
with acoplanarity distributions extracted from in- 
dependent samples of 1 ab~^ luminosity each, 
with a nominal value of (j) = n/A. A precision 
on (f> from such a pseudo-experiment of approxi- 
mately 6° can be anticipated. 

4. Summary 

The combination of generators for production 
and decay of intermediate states, require careful 
treatment of the spin degrees of freedom. In some 
cases one can restrict spin states to pure helicities; 
then generation of intermediate states for individ- 
ual particles can be performed first, and decays 
of each individual particle can be performed later. 
The general case, when full quantum mechanical 
spin correlations are included was also discussed. 
Technical constraints for the solution based on 
kinematical information provided by the produc- 
tion programme to be used by the decay routines, 
were presented. In this context gramatic rules on 
how event records are filled in were discussed as 
well. 

Finally discussion of observable for the Higgs 
boson parity measurement at LC, based on such 
a technical solution was presented in detail, as an 
example. It was shown, on the basis of careful 
Monte Carlo simulation of both theoretical and 
detector effects that with the typical parameters 
of the future detector and Linear Collider set-up 
the hypotesis of the admixture of pseudoscalar 
coupling to the otherwise Standard Model 120 
GeV Higgs boson can be measured up to 6° error 
on the mixing angle. 
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Figure 2. The acoplanarity distribution (angle 
if* ) of the decay products in the rest frame 

of the p^ p^ pair. Gaussian smearing of it's and 
Higgs boson momenta, are included. Only events 
where the .signs of the energy differences yi and y2 
are the same, if calculated using the method de- 
scribed in Ref. TF/ and if calculated with the help 
of the T impact parameter Ref. '20J, are taken. 
The thick line corresponds to a scalar Higgs bo- 
son, the thin line to a mixed one. The left figure 
contains events with yiy2 > 0, the right one is for 
yiy2 < 0. 
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